Scaling in directed dynamical small-world networks with random responses 
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A dynamical model of small-world network, with directed links which describe various 
correlations in social and natural phenomena, is presented. Random response of every site to the 
imput message are introduced to simulate real systems. The interplay of these ingredients results 
in collective dynamical evolution of a spin-like variable S(t) of the whole network. In the present 
model, global average spreading length {L) s and average spreading time (T) s are found to scale 
as p~ a lnN with different exponents. Meanwhile, S behaves in a duple scaling form for N N*: 
S ~ f{p~®q 1 t), where p and q are rewiring and external parameters, a,/3, 7 and f(t) are scaling 
exponents and universal functions, respectively. Possible applications of the model are discussed. 
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Many natural, social and economic systems, are 
well described as networks in which sites represent agents 
and links represent correlations between them (see [1,2] 
for reviews). Among networks with complex topology, 
small-world networks (SWN) [3] are based on a highly 
connected regular lattice in which a fraction p of origi- 
nal links between sites are replaced by randomly rewired 
links. It has played an important role in the study of 
real processes [4] , such as various problems in the held of 
physics [5], and those outside it [6-10]. 

As a platform for understanding processes in sys- 
tems where SWN topology crucially determines their dy- 
namics, Watts [4], Lago-Fernandez et al. [6], Strogatz [9], 
and, Kuperman and Abramson [10] have analyzed several 
types of real problems. Activities of individual agents 
(sites) in their models evolve in time, and couple with 
each other to induce system response from global coher- 
ence, while connections (links) of networks are kept un- 
changed throughout the process. However, interactions 
in most social or biological systems are not static and 
change with time from common knowledge. When states 
of sites at the ends of a link change with time, correlation 
between them could correspondingly change. This effect 
can result in novel dynamical scaling properties which do 
not appear in previous works. 

Suppose an exclusive shop sells a new brand of 
portable commodity without advertising, this message 
can be spread out by showing it in the population, and 
every new buyer can be a source of message, which change 
his/her role in the crowd in this sale process. Mean- 
while, interacting relations among people change contin- 
uously, they are not hxed due to everyday activities. Sim- 
ilar things can occur when a contact-infectious disease 
with all possible incubation (infected-but-not-infective) 
time propagate through an acquaintance network start- 
ing from an index patient. When an infected one finally 
becomes infective, he changes his role to all others. Some- 
times people in incubation periods continuing their social 



activities even are unaware of infection. 

In most information or disease diffusion processes 
effective interactions are asymmetric or even unidirected, 
which has been investigated in a lot of literature [11-14]. 
The direction of a link between a pair of given sites de- 
pends on their states. For example, infection passes from 
infectives to susceptibles, which defines an active link 
from the former to the later. Meanwhile, any link be- 
tween two infectives or two susceptibles is inactive in an 
epidemic spreading. As another example, some special 
cases of food webs which show small- world effect [15,16], 
can be described with inherent directed trophic links [11] 
and their structures can vary with time [12]. These mod- 
els can illustrate how human pollution or disaster on the 
environment to be transported up through species at dif- 
ferent trophic levels. Recently a non-equilibrium phase 
transition is displayed in two-dimensional directed SWNs 
[17]. 

Usually, individuals in real systems react to input 
message not simultaneously [6,10]. Incubation time of 
people for the same infectious disease varies depending 
on their resistivity to it. In parallel, delaying interval 
from getting a message of the commodity to buying it 
also vary according to the spending potential of a par- 
ticular buyer. In one word, there always exist stochastic 
relaxations or random responses in such dynamical pro- 
gresses. 

In this Letter we attempt to point out that three 
ingredients mentioned above — dynamically varied inter- 
actions, directionality of links, and random responses of 
agents can interweave together in a real system, and their 
interplay motivates the present work. We combine these 
ingredients to form a dynamical directed small- world net- 
work (DDSWN) on which we simulate practical spread- 
ing processes. With the model dynamical scaling prop- 
erties of message diffusion and collective response to the 
message are revealed by numerical simulation and data 
analysis, which have not been observed from static undi- 
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rectional SWN models previously. 

The network can be built with the following rules: 
(1) A rewired one-dimensional (ID) circular lattice with 
coordination number 2z is constructed. A site labelled 
j is randomly chosen from all the N sites as the only 
initial seed in spin-up state, while other sites of indices 
i =/= j arc in spin-down state. So the seed is the first one 
that is able to send out the message to other sites and 
to convert spin states of receiving sites. (2) Consider- 
ing the discrepancies among sites in the response to the 
message, we define the response time for every site i as 
the relaxation period from the receiving of the message to 
the flipping into the spin- up state. The Tj's are random 
variables obeying distribution probability P(rj). In this 
paper we consider two types of probability the uniform 
and Poisson distributions. The former can be written as 



P(n) = -0(tO% - 7-0, 



(l) 



where 6(x) is the step function, being for x < and 1 
for x > 0, and q is a parameter characterizing the width 
of the distribution. Here Ti is uniformly distributed be- 
tween and q, so q/2 is the average response period. For 
the latter, we use the discrete Poisson distribution 



P(n) 



q k e 



2 /k\, 



(2) 



where k is the integer part of n, and q corresponds to 
the average value of Tj. In both distributions parame- 
ter q measures the average and diversity of the random 
response periods. (3) We suppose that the spatial con- 
figuration of Ti is time independent, but connections in 
SWN are updated in every time unit To, i.e., a SWN with 
new links is produced at every moment tu — Mt with 
M being a positive integer, from the regular ID circu- 
lar lattice by rewiring p fraction of bonds. So we have 
dynamically varied SWNs with the same rewiring proba- 
bility p. (4) All sites evolve their spin states in a parallel 
way. Once a site of spin-down state is at the end of its 
relaxation time after receiving the message, it flips onto 
spin-up state and keeps it, being capable of sending out 
the same message. Therefore, the links in this dynamical 
SWN are directional as the message can only be trans- 
mitted from a spin-up site to a spin-down site, and the 
transmission from a spin-down site or to a spin-up site 
is ineffective. Actually, the direction of a link reflects 
the irreversible nature in diffusion processes, such as the 
spreading of pollution in food webs [9,11] and diffusion 
in directed percolation models [18]. 

For an undirected and static SWN, the small- 
world effect can be described by the characteristic path 
length I, i.e., the averaged shortest distance between any 
two sites. It follows the scaling law [1,2], 

l(N,p) ~ (N*y/ d F(N/N*), (3) 

where the scaling function F(u) has the limits F(u) ~ 
u 1 /" 1 for u < 1, F(u) ~ lnu for u > 1, and N* ~ p^ 1 



is a crossover size separating the big- and small-world 
regimes. In DDSWN, however, not all sites are effec- 
tively connected with each other since only links from 
sites in spin-up state to those in spin-down state are 
active. Therefore, I loses its original meaning in the 
present model. In order to investigate how wide and 
how fast the message spreads out in DDSWNs, we in- 
troduce the time-dependent average spreading length, 
(L/j) = Ylv=i n vl v / Z)£=i n v, and the average spreading 
time, (Tp) = YTv=\ n v^v I YTv=\ n v ^ where \i is the time 
index, t M = (ito, is the number of receivers created 
between the (/z — l)th and the /zth moments, and is 
the average length from the initial seed to these receivers. 
We calculate them by numerical simulations in which the 
data are obtained by ensemble average over at least 50 
random configurations with different {t^} and location of 
the initial seed. The obtained results are shown in Figs. 
1 and. 2. 
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FIG. 1. Log-log plot of saturated average path length (L) s , 
in the case of uniform distribution of n, versus p. Two lines 
have the same average slope, 0.86. Upper-left inset: (L) as 
a function of t for the curve with symbols (+) in the main 
panel. Lower-right inset: log-linear relation of (L) a versus N. 
Other parameters for these curves are z = 5 and q = 10, and 
the.tti mc is in units to. — 1 1 




FIG. 2. The same as those in Fig. 1, but for Poisson distri- 
bution of Ti . In this case the average slope of lines in the main 
panel is 0.84. q = 10 and other parameters are indicated in 
the figure. 
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It is found that for given N and z the average path 
length (L(t)) (here we indicate the time-dependence ex- 
plicitly and omit the suffix p) increases monotonically in 
time towards its p-dependent saturation value (L) s (see 
the upper-left insets of Figs. 1 and 2), which shows the 
finite-size effect of p. Summation over n v in the calcula- 
tion of (L) a tells us that eventually all sites will be suc- 
cessively visited by the message sooner or later, suggest- 
ing the higher adaptability of the dynamical models in 
description of the diffusion processes compared with the 
static directed SWN models [14]. We will use (L) s and 
(T) s to characterize global properties of the networks, 
since in saturation the quantities are no longer time de- 
pendent. Compared with the scaling relation Eq. (3) for 
l(N,p) in the static SWNs, the corresponding quantity 
(L) s in the DDSWNs scales with \nN for given p (see 
the lower-right insets of Figs. 1 and 2), but it depends 
on p in a power law (L) s ~ p~ aL for a fixed N (see the 
main panels of Figs. 1 and 2, where the slope of straight 
lines gives the value 1 — at). Therefore, we have a new 
scaling formula 

(L> s ~f>- a MnAT, (4) 

where = 0.14 ± 0.02 for uniform distribution and 
Qi = 0.160 ± 0.005 for Poisson distribution, extracted 
from the data in Figs. 1 and 2, respectively. Meanwhile, 
the average spreading time (T) s behaves in the same 
form as (L) s , but with different parameters. We obtain 
a T = 0.29 ± 0.01 and a T = 0.245 ± 0.006 for uniform 
and Poisson distributions, as extracted from the data in 
Figs. 3 and 4, respectively. The message from the initial 
seed passes through tree-like paths, and (L) s and {T) s 
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FIG. 3. Saturated average spreading time (T) s as a func- 
tion of p for the uniform distribution of the response time. 
Other parameters are the same as those in Fig. 1. The aver- 
age slope of the lines in the main panel is 0.81. 

are the maximal path length and spreading time by which 
the message reaches the final sites. The scaling relations 



of (L) s and (T) s are unique features of DDSWN model. 
The simulations suggest that the p-dependence of the 
saturated (L) and (T) is attributed to the dynamically 
varied connections and the finite-size effect, the factor 
In N in the scaling relation can be attributed to the geo- 
metrical characteristics of the SWN background, and the 
factor p~ aL is originated from the diffusion mechanism of 
the directional links, distinctive in comparison with the 
behavior of / in the ordinary SWNs. The exponents oil 
and cut still depend on specific form of the distributions 
of the relaxation time. On the other hand, as it should 
be anticipated, (L) s ~ N for N <C N* in the big- world 
regime. 
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FIG. 4. The same as those in Fig. 3, but for the Poisson 
distribution of the response time. Other parameters are the 
same as those in Fig. 2. The average slope of the lines in 
main panel is 0.755. 

We now turn to the behavior of S(t), defined as the 
number of all sites in spin-up states at moment t divided 
by N. It represents the ratio of the sites in spin-up state, 
and its temporal derivation, S'(t), reflects the probability 
of a spin-down site to be flipped into the spin-up state 
in a unit of time at moment t. Because in this model 
all the sites should be in spin-up state after an enough 
time, S(t) — ► 1 at the limit t — > oo. We numerically cal- 
culate S(t) for various values of parameters. The results 
for different values of p and fixed q in the case of uniform 
distribution are plotted in Fig. 5(a). All the curves can 
merge together if the time is rescaled as t sc = (p/po^t, as 
can be seen from Fig. 5(b) in which we use po = 0.1 and 
[3 = 0.20 ± 0.005 is extracted. The time rescaling is done 
for the Poisson distribution, as shown in Figs. 5(c) and 
5(d), and we get (3 = 0.186 ± 0.006. At the same time, 
if q is varied as well, we can also use the time rescaling, 
^sc = (?/9o) ~ 7 i sc to merge the curves, as shown in Fig. 
6. From this we extract the exponent 7 = 0.750 ± 0.003 
for the uniform distribution and 7 = 0.880 ± 0.006 for 
the Poisson distribution, both with q = 10. Thus, we 
have the following scaling form of S(t): 
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S(t) ~ S(p~^t) EE f(t), 

where t = Po~^go 7 4c an d /(*) ^ a universal function for 
a given distribution. 
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(5) example, usually have a few directed links [16] from the 
basal species to the top ones, and moreover, have average 
trophic distance scaling in In TV [19]. All of them share 
typical values of (L) s in the DDSWN model. This im- 
plies that the model can capture the essential feature of 
dynamical directionality of food webs although the small- 
world effect in this problem is still in dispute [15,16]. By 
the way, on revision of our submitted manuscript, simi- 
lar In AT scaling behaviors of SWNs were independently 
found in quantum and classical diffusion problems [20]. 

This work is supported by China Aeronautic Sci- 
ence Fundation under Grant No.00J52079. 



FIG. 5. The ratio of spin-up sites versus time t (a and c) 
and versus rescaled time t ac = (p/po) 8 t (b and d) for the uni- 
form (a and b) and Poisson (c and d) distributions, po —0.1 
is used. Other parameters are: N — 2000, z = 5, and q — 10. 





FIG. 6. The ratio of spin-up sites versus t sc (a and c) and 
versus t' sc = (q/qo)~''t sc (b and d) for Poisson (a and b) and 
the uniform (c and d) distributions, go = 10 is used. Other 
parameters are: N = 2000 and 2 = 5. 

The distributions of the response time used here 
are phenomenological functions reflecting the nature of 
the DDSWN. Other forms can also be used to describe 
specific systems. The results of our simulations with 
different distributions suggest that there exists common 
scaling behavior of the diffusion dynamics in DDSWNs 
and only the values of exponents depend on the form of 
distribution. For a specific process, the response time 
has a definite meaning. Sales problem of a new type of 
commodity, for example, r, reflects the delaying time be- 
tween the awareness and the buying of the commodity for 
a consumer, and q is related to the discrepancy in con- 
sumptive potential among them. Food webs, as another 
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